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Abstract 

We consider a dilute homogeneous mixture of bosons and spin-polarized 
fermions at zero temperature. We first construct the formal scheme for car- 
rying out systematic perturbation theory in terms of single particle Green's 
functions. We especially focus on the description of the boson-fermion in- 
teraction. To do so we need to introduce a new relevant object, the renor- 
malized boson-fermion T-matrix which we determine to second order in the 
boson-fermion s-wave scattering length. We also discuss how to incorpo- 
rate the usual boson-boson T-matrix in mean-field approximation to ob- 
tain the total ground state properties of the system. The next order term 
beyond mean-field stems from the boson-fermion interaction and is pro- 
portional to obf^- The total ground-state energy-density reads E/V = 
+ cb + {^-nTl^ a-QY^Bn-p / m)\i + aBF^F/('5)/7r]- The first term is the kinetic 
energy of the free fermions, the second term is the boson-boson mean-field 
interaction, the pre-factor to the additional term is the usual mean-field con- 
tribution to the boson-fermion interaction energy, and the second term in the 
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square brackets is the second-order correction, where f{S) is a known func- 
tion of (5 = (mB — mF)/(mB + mp). We also compute the bosonic and the 
fermionic chemical potentials, the compressibilities, and the modification to 
the induced fermion-fermion interaction. We discuss the behavior of the total 
ground-state energy and the importance of the beyond mean-field correction 
for various parameter regimes, in particular considering mixtures of ^Li and 
^Li and of ^He and ^He. Moreover we determine the modification of the in- 
duced fermion-fermion interaction due to the beyond mean-field effects. We 
show that there is no effect on the depletion of the Bose condensate to first 

order in the boson- fermion scattering length obf- 
PACS numbers: 03.75.Fi, 03.70.+k , 01.55.+b 
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I. INTRODUCTION 



Following the spectacular success in achieving Bose-Einstein condensation in trapped, 
dilute atomic gases in 1995 Jl]-^, there has been an explosion of experimental and theoretical 
activity on this newly accessible state of matter (for recent reviews focusing on different 
experimental and theoretical aspects see for instance |0-§[). More recently, there has been 
increasing interest and experimental activity also in quantum degenerate ultra-cold Fermi 
gases [P|-P^, in particular because of the possibility of observing a BCS type transition 
in a dilute atomic gas [^,0 • Dilute mixtures of ultra-cold gases of bosonic and fermionic 
atoms are also receiving increased attention, in particular because sympathetic cooling of the 



fermions by the bosons is an important means of their achieving quantum degeneracy |Tl|-[14 



and also because bosons can mediate an induced (attractive) fermion-fermion interaction 
T7| . Moreover mixtures of atomic ^He and ^He have become interesting in their own right 



after the recent achievement of Bose-Einstein condensation in metastable ^He |jT8|,|T9[, as 
they could represent a bridge towards the understanding of superfluidity in helium. 

Current analyses of dilute mixtures of ultra-cold atomic boson and fermion vapors are 
based on mean-field approximations. They include for example the work on stability con- 
siderations for homogeneous systems by Viverit et al. |T^, and the calculation of density 
distributions and phase separation of trapped mixtures by Nygaard and M0lmer . Some 
interesting effects have been studied by Bijlsma et al. where an effective modification 
of the fermion-fermion scattering length, mediated by boson-fermion scattering processes, 
was determined. Pu et al. determined the phonon spectrum of the Bose condensate in a 
boson-fermion mixture at zero temperature ||22|| . 

Mean-field approaches have proved to be extremely useful in the theoretical and ex- 
perimental study of Bose-Einstein condensed dilute atomic gases, and are likely to prove 
similarly useful for quantum degenerate mixed boson-fermion systems. It is nevertheless de- 
sirable to consider beyond mean-field effects, and under what circumstances they are likely 
to be most relevant. For pure (unpolarized) fermion p3|-p^ and pure boson [pl|-pl[] systems. 
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expansions of the ground-state energy, in terms of the small parameters kpapp and y t^bObb 
{kp is the Fermi wavenumber, ne the boson density, app and obb the fermion-fermion and 
boson-boson scattering lengths), are well established. These expansions go beyond mean- 
field approximations while still depending only on the s-wave scattering lengths. Although 
determined for homogeneous systems, the use of beyond mean-field corrections arising from 
consideration of such expansions may be readily extended to the experimentally relevant 
case of inhomogeneous trapped gases by application of the local density approximation. In 
general, the beyond mean-field corrections for the bosons are smaller than for the fermions, 
since the exponent of the small dimensionless parameter n^^^a (n is the density parameter 
and a the scattering length) is 1 in the fermion case but 3/2 in the boson case. 

In the case of dilute fermions immersed in a Bose gas an expansion of the ground-state 



energy in terms of the small parameters J cibb^b and np/n-B, where rip is the fermion density. 



was performed by Saam [32]. This was motivated by considering quantum-degenerate dilute 
gases as a model for the behavior of superfiuid helium, where the assumption of n-p/n-Q as a 
small parameter is justified by the much greater natural occurrence of bosonic ^He compared 
to that of the fermionic ^He isotope. 

Systems where there are vastly more bosons than fermions are certainly experimentally 
achievable in dilute atomic gases, and it can in fact be advantageous to have an excess of 
bosons in order to enhance sympathetic cooling [|lT|. However, there is in principle no a 
priori reason to confine theoretical analyses to such systems. In fact, in recent experiments 
13| , p^ the numbers of fermions and bosons are comparable. Thus motivated, in the present 



paper we derive a systematic perturbative expansion for the ground-state energy and other 
related relevant physical quantities for dilute Bose-Fermi mixtures at zero temperature and 
for arbitrary ratios of the boson and fermion densities. In this way we determine the lowest- 
order correction to mean-field in the case of weakly interacting bosons and spin-polarized 
fermions in terms of the the gas parameter kpa-^p, where obf is the boson- fermion s-wave 
scattering length. The ground-state energy thus derived can then be implemented, in local 
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density approximation, as the energy functional for the study of the experimentally relevant 
case of trapped mixtures, in complete analogy with the pure bosonic and pure fermionic 
cases. 

The plan of the paper is as follows. In Section y we introduce the basic Hamiltonian for 
a system of interacting bosons and spin-polarized fermions, expressed in its grand-canonical 
form after performing the Bogoliubov replacement. In Section |T|we define the one particle 
Green's functions needed for a systematic field-theoretical analysis of the boson-boson and 
boson-fermion interactions, and we determine the associated Feynman rules. In Section 



rVl we implement the perturbative expansion by introducing the boson-fermion self-energy 
and the renormalized boson-fermion T-matrix in ladder approximation and by solving the 
corresponding Bethe-Salpeter equation to second order in /cf^bf- In Section we exploit the 
results obtained in the previous Sections to compute some relevant physical quantities. In 
particular we provide the expression for the ground-state energy density to second order in 
the gas parameter, the bosonic and fermionic chemical potentials, the compressibilities, and 
the induced fermion-fermion interaction. We then compare the results thus obtained with 
actual and foreseeable experimental situations to assess the relative importance of higher- 
order corrections with respect to the mean-field results. In Section [V^ conclusions are drawn 
and some possible future developments are discussed. 



II. SYSTEM 
A. Hamiltonian and ground-state energy 

1. Many-body Hamiltonian 

We consider a homogeneous mixture of interacting bosons and fermions, imposing peri- 
odic boundary conditions on a volume V. In complete generality there are thus boson-boson, 
boson-fermion, and fermion-fermion interactions to consider. However, for spin-polarized 
fermions, there is no s-wave scattering contribution to the fermion-fermion interaction . 
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The first non-vanishing contribution is due to p-wave scattering, which can generally be 
neglected when compared to the boson-boson and boson-fermion interactions, which are 
due to s-wave scattering. We thus take into account s-wave scattering between bosons, and 
between bosons and fermions only. 

In second-quantized form, the Hamiltonian describing this situation is 



H = fB + fF + U + V, (1) 



where 



ft2 ^ 

Tb = / rf3xV$t(x) • V$(x), (2) 

Tf = / d^xV-^^x.) ■ V^(x), (3) 

U = J J rfWx'<|t(x)^t(x')f/(|x - x'|)^(x')l>(x), (4) 
V = lj J d3xd3x'<l>t(x)<|t(x')y(|x - x'|)<l>(x')<l>(x), (5) 

and where $(x) is a bosonic field operator, \l/(x) is a fermionic field operator, and me and 
mp are the respective masses of the bosons and fermions. For later reference we also define 

Ho = fB + fp, (6) 
W = U + V. (7) 



2. Mean field theory 



It is straightforward to determine a zero-temperature mean field theory for H . Em- 
ploying the well-known Thomas-Fermi approximation, the mean field ground-state energy 
density is 

E 3 h^kp 27r^^aBF 27r^^aBB o / \ 

— = -7^ riF + ubUf + n^, (8) 

V 5 Zmp m rriB 

where m = mYm-Q/{m-p + itl-q) is the reduced mass, and k-p = {Q'n'^n-pY/'^ is the fermi wave- 



number P^. In the case of a pure fermionic system, beyond mean-field corrections to the 



ground state energy-density are given by [p3 



5mF 



Up 



. 128, ,2 
1 + -—kpap + [kpaFj + 



(9) 



For pure bosons corrections to the ground state have been calculated by, e.g., Hugenholtz 
and Pines ^8 and by Wu [29|. These corrections are obtained via a perturbative expansion 



in terms of the bosonic gas parameter y^t-b^Ib- already mentioned, this parameter is 
in general smaller than the fermionic gas parameter (see also Sec. V B ). Our goal is thus 
to determine a general expression equivalent to Eq. (^, taking into account boson-fermion 
interactions, while neglecting corrections proportional to higher powers of the bosonic gas 
parameter. 



B. Bogoliubov replacement and grand-canonical Hamiltonian 

In order to determine the energy functional to higher order than in Eq. (H), we will adopt 
a perturbative approach using one-particle Green's functions, in a way essentially equivalent 
to the field-theoretical treatment of pure bosonic and fermionic systems [^,^]. We thus 



first carry out the Bogoliubov replacement [^], where the condensate bosons are treated as 
a c-number field: 

$(x) = v/^ + 0(x), (10) 

where uq = Nq/V is the condensate density, and Nq is the number of (condensate) atoms 
in the A; = mode. This prescription breaks particle number conservation (see PB|-PU| for 
alternative Bogoliubov replacements that preserve particle number conservation); average 
particle number conservation is assured by introducing the grand-canonical Hamiltonian 

k = H-flBNB, (11) 

where /xb is a Lagrange multiplier, to be identified with the boson chemical potential p5 . 
Substituting Eq. (p!OD into Eq. (|ri|), the grand-canonical Hamiltonian reads: 

K = Ko- I^bNo + Ui + U2 + U3 + V, + V2 + V^ + V^ + V5 + Ve + Vj, (12) 



where 



Ko^^J c^3xV0t(x) • V0(x) + ^ / c^'xV*t(x) • V*(x) - /xb / dW(x)0(x), (13) 

U, = noJ I d'Wx'^t(x')[/(|x - x'l)^(x'), (14) 

C/2 = ^ y y'ci3^(iV*t(x')C/(|x-x'|)*(x')0(x) +h.c., (15) 

Us^ J J c^WxVt(x)^t(x')[/(|x - x'|)^(x')(^(x), (16) 

14 = ^^0 / / ^'xdVV^dx - x'l), (17) 

% = noV^ J J d=^x(^VT/(|x - x'|)0(x) + h.c, (18) 

^no J J ci^x(iV\/(|x - x'|)(^(x')(^(x) + h.c, (19) 



V4 = noj j c^Wx>t(x)\/(|x - X |)</)(x'), (20) 

T/5 = no/ y'ci3xciV0t(x)v^(|x-x'|)0(x), (21) 

V^6 = v^/ / rfWx>t(x)\/(|x - x'|)0(x')0(x) + h.c, (22) 

Vr^lj J ciWx'0t(x)0t(^')^(|x _ x'|)0(x')0(x). (23) 



III. SYSTEMATIC PERTURBATION THEORY WITH GREEN'S FUNCTIONS 

A. Green's functions: definitions 

The boson (B) and fermion (F) Green's functions for the boson-fermion system are 
defined as 

zGB(x,t,x',0 = (e|T[|.(x,t)|.t(x',t')]|0, (24) 
iGp(x,f,x',i') = (e|r[*(x,i)*t(x',i')]|e>, (25) 

where the time argument in $(x, t) and ^(x, i) means they evolve according to Heisenberg's 
equations of motion, T denotes the time ordered product, and |^) is the ground state of K 
(we similarly define |^o) to be the ground state of Kq). We use the Bogoliubov replacement 
to write 
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zGbIx, t,x',t') = no + iG'^{x,t,x',t'), (26) 

where 

tG'^{^,t,^',t') = {m4>M4>\^',t')]\0 (27) 

is the propagator for the non-condensate bosons. 

B. Perturbative expansion 



The Green's functions can be evaluated in perturbation theory where W is the 
perturbation to Kq. Thus 



where 



zGirHx,t,x',t') = (eo|T[5(")0(x,t)0t(x',t')]|^o), (30) 
zGirHx,t,x',t') = (eo|T[5(")§(x,t)^t(x'^i')]|^^)^ (31) 
5^"^ = ^ (x) V • • • / ^'^nT[H^(ti) . . . W{Q]. (32) 



Operators with a tilde are defined to be in the interaction picture, i.e. 0{t) = 
exp{iKot/h)0 exp{—iKQt/h). In the limit of a non-interacting system {W 0) the Green's 
functions reduce to the zeroth order terms in the expansions, so that 

(x,t,x',t') = z4°)(x,t,x',t') = {^o\T[4>{^,t)4>K^',t'mo), (33) 
zG°(x,t,x',t') = zGf)(x,t,x',t') = {^o\T[^i^,m\^',t'mo). (34) 



C. Evaluation of terms using Wick's theorem 



Equations (|5DD, (pi|), and ( P^D can be evaluated by Wick's theorem, which states that 
the vacuum (non-interacting ground-state) expectation values of time ordered products of 
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operators can be expressed as the sum of all products of contractions of pairs of operators 
in the time-ordered product The contraction of two operators is defined as 

d{tf^P{t'f^ = T[d{t)P{t')] - : d{t)P{t') : , (35) 

where : 0{t)P{t') : is the normal ordered product. In particular, 

0(x,t)«0t(x'^i')« = 4>\^',t'Y^4>(^,tf) = zG°(x,t,x',t'), (36) 
^(x, t)«§t(x', = -§t(x', t')W^(x, tf^ = «G^(x, t, x', t'), (37) 

and all other contractions of pairs of operators G {(/)(x, t), (/)^(x', t'), \l/(x", t"), \&^(x"'', t'")} 



vanish (see also Appendix Substituting Eqs. (|36D and (^) into Eqs. 
fB^), the first order terms can be determined to be: 



and 



—I 



X '-noiG%{yly',)iGl{x^,yn 
-iGl{y^,,y^)iGl{x^,x^,)iGl{xly^^ 
+V{x^, - yt) 

+n,iGl{x^,x'i)iGl{x^^,y^) 
+r^o^G°(x^a;^')^G°(x^x^) 

^[-n,iGl{yly^,)iGl{x^,y^) 
-n,iGl{x^,y'i)iGl{yly^)] 
+V{x^,-y^,)'^iG^{x^,yn] 
J d'x^,d^y^[u{x1 - y^) 



(38) 



n 



(39) 



(40) 



where we have used the more compact four- vector notation [x^ = {t, x)], and defined f/(x^ — 



t/(x-y)5(x°-i/°) and V{x^'-y^') = V{x-y)6{x^ -y^). Note that Gl^(t,x,t,y) = 



[i.e. there are no boson loops at zero temperature) 
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Higher order terms may be similarly evaluated, and will similarly be expressed in terms 
of integrals over products of noninteracting Green's functions, condensate factors no, and 
interaction terms. We represent these graphically (see Fig. |I]): straight lines for fermions, 
wiggly lines for non-condensate bosons, dashed lines for condensate bosons, and zigzag lines 
for interaction terms (whether it is a boson-boson or boson-fermion interaction is clearly 
determined by the kinds of particle lines attached to the vertices of the interaction line): 
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FIGURES 

•^/VWWV* = Uix" - y^) or V{x^' - y^) 
• • = v^- 



FIG. 1. Definition of the diagram lines 



As is usual [^,^,^ , all disconnected graphs in the numerator can be factorized out by 



oo 



the denominator, so that 

G'^{X^, yn = £ 4" l/'^)connected , (41) 
n=0 

oo 

G'F(X^ yn = E l/'^)connccted • (42) 

n=0 

Noting that each connected graph essentially appears n! times, with simple permutations on 
the labeling, when composing such graphs we integrate over all internal variables and affix 
a factor of {i/h)"'{—l)^{—i)'-^, where n is the number of interaction lines, F is the number 
of closed fermion loops, and C is the number of dashed boson lines. 

D. Feynman rules 

For homogeneous systems it is convenient to Fourier transform to energy-momentum 
space, so that: 

^f(p') = n 210 ]r- 1 TV' (44) 

— hp'^/Zmp + tsgn[p — k-Fju 

where sgn(A;) = 1 for A; > and = —1 for A; < (we write p for |p|). The appropriate 
Feynman rules for the boson (fermion) Green's function in this representation are then: 
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1) Draw all topologically distinct connected diagrams with one outgoing external wiggly 
boson (fermion) line and one incoming external wiggly boson (fermion) line, no external 
fermion (boson) lines and no internal dashed boson lines, n zigzag interaction lines, each of 
which is attached at one vertex to an incoming and an outgoing boson line (either wiggly 
or dashed), and at the other vertex either to an incoming and an outgoing boson line, or to 
an incoming and an outgoing (not necessarily distinct) fermion line. Each vertex must be 
attached to exactly one zigzag interaction line. 

2) All wiggly boson lines must run into the same direction and there are no closed boson 
loops. 

3) Each dashed boson line corresponds to a factor of i/no, each wiggly boson line to a 
factor of G^{k^), each fermion line to a factor of G^{k^), each boson-fermion interaction 
line to a factor of U{k^) = U(k), and each boson-boson interaction line to a factor of 
Vik") = V{k). 

4) Assign a direction to each interaction line; associate a directed four-momentum with 
each line and conserve four-momentum at each vertex. Each dashed boson line carries 
four- momentum and each wiggly boson line has four- momentum ^ 0. 

5) Integrate over the n independent four-momenta. 

6) Affix a factor of (i//i)"(27r)~^*^"^(— 1)^(— i)*^, where F is the number of closed fermion 
loops and C is the number of dashed boson lines. 



IV. DETERMINATION OF THE BOSON-FERMION T-MATRIX AND 
SELF-ENERGIES IN LADDER APPROXIMATION 

A. The Hugenholtz-Pines theorem 



According to the Hugenholtz-Pines theorem p8|j4^ , the bosonic chemical potential jiB, 
defined as 
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dE/V 
drin 



(45) 



is given by 



/XB = /iSB(o)-;iSi2(o), 



(46) 



where Eb(0) and ^12(0) are the proper self-energies for the bosons due to their interaction 
with both bosons and fermions, evaluated at p'^ = (in what follows we call them the 
bosonic self-energies). The self-energies are in general related to the Green's functions by 
the Dyson equations. The Dyson equation for the bosons is given by: 



I 



Gim 
Gl{-v^) 



I 



+ 



Gl{v^) 
1^ Gl{-v' 



\ i 



I 



\ 



I 



G2i(p^) G'b(-p^) 



(47) 



Sb(p^) Si2(p'^ 
S2l(p'') Sb(-p'^) 

where have introduced the anomalous boson Green's functions 6*12 (p^) and G2i{p^) (de- 
fined as the Fourier transforms of Gi2{x^,y'^) = (^|7'[0(a;^)0(y'')]|O G2i{x'^,y^) = 
(G|[0"'^(x'^)0"'^(?/'^)]|^), respectively). The Dyson equation for the fermions takes the much 
simpler scalar form: 



GFiP") = GUpn + GUpn^FiP^GFiP' 



(48) 



where Sf(p^) is the proper self-energy for the fermions due to the interaction with the bosons 
(the fermionic self-energy). 

B. The self-energies in ladder approximation 



As we are considering a dilute system, in terms of Feynman diagrams only diagrams 
with interaction lines between two systems of connected propagators are important P^J2H 
(ladder approximation). This is expressed in terms of the boson-fermion and boson-boson T- 
matrices in Fig. p|, where the boson-fermion T-matrix Tbf in ladder approximation is defined 
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in Fig. 1^, the boson-boson T-matrix Tbb (also in ladder approximation) is well known from 
studies of dilute pure Bose systems, and the normal (diagonal) bosonic proper self-energy is 
given by 



Sb(p^) = Sbf(p^) + Sbb(p^). 



(49) 



n^BFip") 



P" A 



p" A 
p" A 



p" A 



p" A 



p" A -p" A 



?!Sj,(p") 



FIG. 2. The self-energies in ladder approximation, expressed in terms of the T-matrices. 



The proper self-energies can thus be determined by adding the proper self-energies of a 
system of interacting bosons to those of a hypothetical mixed boson-fermion system where 
there are boson-fermion interactions only |Q. This result arises from our use of the ladder 
approximation, and is not in general true (there also exist, for example, inseparable three- 
legged "ladders" consisting of a boson-boson and a boson-fermion ladder joined by a common 
boson leg, but these clearly involve three-particle processes). For such a hypothetical mixed 
system, the only self-energies we need to consider and to evaluate are Sbf(p^) and Sf(p'^), 
which can be written algebraically as: 
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Pi 



P3 . 



IP5 



(27r)4 



dVTBF(p^F,p^F)G'°(F), 
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(50) 
(51) 



FIG. 3. The boson- fermion T-matrix. 



C. Bethe-Salpeter equation for Tbf 



Pi a P2 




p5*A Ap^ 
•WVWtAM* 

P^iip'^-Pl kpl 




FIG. 4. The integral equation for Tbf- 



The boson-fermion T-matrix Tbf can also be represented recursively, as shown in Fig. ^. 
If we now transform to center-of-mass coordinates, 
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(52) 



the algebraic form of the equation represented in Fig. reads: 



T^^ikl fc^, Pn = f/(ki - k2) + j-^^ J d'kUik, - k) 



(53) 



This is a kind of Bethe-Salpeter integral equation, which we will now solve recursively for 
low momenta, stopping at order a^p . As the interactions are instantaneous, the only 
frequency dependence in T^piki, kl^, P^) is in P° [0,^. Thus, a contour integration over 
fco in Eq. (m yields: 



TBF(ki,k2,P^) = f/(ki-k2 



(2vr) 



Uiki - k)TBF(k, k2, P^)9{\P/2 - k| - kp) 



(54) 



"/iPO - n'^(P/2 + k)2/2mB - n'^(P/2 - k)V2mF + + iu 
We now express Eq. (|5^ in terms of the free scattering amplitude /(ki, k2), first by formally 



inversion (see |p4|| ): 
27ih^ 



m 



/(ki,k2) = f/(k2-ki 



1 



^,^Uik2-k)2nh'f{k^,k)/m 



(27r)3y ^ " h\j/2m - h\y2m + iu ' 
and then by exploiting the resulting expression to rewrite Eq. (0) as: 

rBF(ki, k2, Pn = ^/(k2, ki) + / d'k^fik, kOTBF(k, k2, P^ 



(55) 



X 



^(|P/2-k| - A;f) 



hPO - r(P/2 + k)V2mB - r(P/2 - k)V2mF + + 
1 



(56) 



/i^kf /2m - ;i^kV2m + iz/_ 

For low momenta the vacuum scattering amplitude /(ki,k2) can be expanded to second 
order in the scattering length Obf (see [P5[): 

/(ki, k2) ^ obf - ialpk, (57) 

where k = ki = k2 ^ 0. We insert this into Eq. (|56D , iteratively substituting Eq. (|56D into 
itself, and consistently keeping terms up to quadratic order in obf only. This produces 
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TBF(ki, k2, P^) ^ —[aBF - ^alph] + 



m 

X / d^k 



e{\p/2-k\ - kp) 



- h^{P/2 + k)V2mB - n^{P/2 - k)V2mF + /i + zi/ 
1 



(58) 



h^'kl/2m-h^k^/2m + iv_ 

the renormalized second order expansion of the boson-fermion T-matrix. The integral can 
be evaluated (see App. 0) to give 



TBF(ki,k2,P^ 



27111^ 2h^alpkF 



m 



-a-BF + 



m 



a'i-ph'^ f m^kl m^P ^ r— me-DA , k^ + mP Im-Q + \fD + ii' /2a-QF\fD 
' 2my D m 



2m2 \ P 



2m? 



where 



me 



P J k^ — mP/m-Q — \/15 ~ iv /2a-QY^fD 



P 



m P I — mB-D\ , k-p — mP m-B 

h 2mv D m 

•rriB P kY + mP/mB 



D + iu/2aBFVD 



y/D-iu/2aBFVD 



D 



m 



niB + m-p 



-P2 + 



2mP^ 2m/i 



h 



+ 



(59) 



(60) 



V. PHYSICAL QUANTITIES 



A. Bosonic chemical potential 



Substituting Eq. (|^) into Eq. (|50|) the equation for Sbf(p'^) can be rewritten as 
- /.4„. ?^BF[(p-q)/2,(p-q)/2,p/- + g^' 1 



{2%)^ ~ ha^ /2m-p + zsgn(g — kY)v 



(61) 



To evaluate this, we substitute Eq. (^) into Eq. (|6TD, and first carry out the frequency 
integral. As the pole in the complex g°-plane of the integrand in Eq. (^) is below the real 
axis, in order to get a non vanishing result the pole of — hq'^/2'mp — isgn{q — kp)u]~^ must 
be above the real axis {q < kp). The frequency integral is thus readily solved by contour 
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integration. The k integration in (|5^) is then very similar to that leading to Eq. (|59D. The 
resulting expression for hJ^-QFip'^) is then: 



(27r) 



p-q p-q 



P 



p + q 



(62) 
all terms 



2 ' 2 ' V 2mF 
We wish to similarly solve this integral to second order in obf- In Eq 
which depend on D have a pre-factor a^p. Thus, in order to get a result for Eq. (p2| ) that 
is correct to second order in obf, it is sufficient to use the zeroth order expression for D. 
Specializing to the case where p'^ = this can be written as 



ml 



I , ^2-9^• (63) 

We now substitute for D in Eq. (|59D, and, after a straightforward (if lengthy) integration 
over q, arrive at 

flBF^F 



2T:fi?aF,F 

^Ebf(O) = n-p 



m 



1 + 



TT 



-m 



where 



(64) 



(65) 



45 ' 8(52 "^1-5' 

S = (mB—m-F)/(mB+mp), and we have used (27?)""^/ d'^k'^GF(kf^) = Gf(x'^,x'^) = inp- Note 
that in this integration we need only consider the real part of the boson-fermion T-matrix, 
as within the range of the integration the imaginary part is zero (see App. 0). The necessary 
expression for the T-matrix is then just given by Eq. (p9|), where we take the absolute values 
of the arguments of the logarithms and set = 0. From the Hugenholtz-Pines theorem [Eq. 
(H)]: 

/xb = ^Sbf(0) + ^Sbb(0)-/iSi2(0). (66) 

Thus, using the expression for I1bf(0) in Eq. (0), and the results from for Sbb(O) and 
Si2(0) (neglecting corrections of the order of the boson gas parameter), 

(67) 



Aiih obb 2-7?^ obf 
/iB = n-B H uf 



1 + ^/(5) 



7r 



rriF, m 

This is exactly equivalent to adding /iSbf(O) to the standard mean-field result for the bosonic 
chemical potential for a pure, self-interacting bosonic system. 
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B. Ground state energy density 



To obtain the ground state energy we simply integrate Eq. (^5]): 



E 
V 



fi{nB)dnB + C(nF), 



(68) 



where C(nF) is a quantity that can depend on the fermion density np only. Considering the 
limit Obf —>■ 0, we see that C{n-F) can only be the kinetic energy for free fermions (the Fermi 
energy density ep) that is: 



CM = 6F = 

5 2mp 



(69) 



Substituting this and Eq. (|67| ) into Eq. (|68D , and then integrating, gives, finally: 



— = ep + cb H npUB 

V m 



1 + 



TT 



(70) 



where eB = 27r/i^aBB'^B/'^B is the boson mean-field energy density. Eq. (|70D is the main 
result of this paper, being the desired extension of the mean-field result Eq. (H). 
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FIG. 5. Plot of f{S), where 6 = {uib — mp)/{mB + rnp), proportional to the correction to 
second order in a = /cfobf to the energy-density functional [Eq. (74)]. The relevant values of f{S) 
for mixtures of ^He and ^H, ^Li and ^Li, ^He and '^He, and '^''K and ^^Rb are indicated. 

It is illuminating to describe Eq. ([7D| ) in terms of the dimensionless gas parameters and 
the dimensionless ratio of the boson and fermion densities: 
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a = a-QF^F, (71) 



(3 = ^nBa^B, (72) 
V = —, (73) 



so that 

E 



2/3 



1 + -fis) 

TV 



(74) 



The corrective term to second order in a is proportional to the rather comphcated function 
f{S), defined in Eq. (|65D, of the relative mass ratio 6; the value of this function will thus vary 
considerably depending on the masses of the atomic species used in any given experiment. 
In Fig. ^ the values for mixtures of ®Li and ''Li, and and ^^Rb, corresponding to real 
experimental configurations currently under investigation, are plotted, as well as that for a 
mixture of '^He and ^He, also a likely candidate for future investigation in ultra-cold dilute 
gas experiments. The value for a hypothetical mixture of ^He and is also shown, as it 
is almost exactly the maximum possible. The function f{6) is always positive in the total 
range [—1, 1] of variations of 6. Note that in the limit m-Q/mF —* oo, one has 6^1 and 
f{6) —>■ 0. Thus the second-order correction to the boson-fermion interaction energy and the 
total boson-boson interaction energy disappear. This is because if the bosons are infinitely 
massive (compared to a fixed, finite fermion mass), then it is impossible for them to be 
scattered out of the condensate, and only the boson-fermion mean field interaction remains, 
since all the bosons can be treated as the (condensate) mean field. In the opposite limit of 
m^/mp 0, the situation is different, because of the Pauli exclusion principle. 
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FIG. 6. Boson-boson (solid line) and boson- fermion (dotted line) mean field contributions 
(a) and second-order corrections (b) to the ground-state energy-density of a ^Li, ^Li mixture, in 

1/3 

units of ep, to the ground-state energy-density in units of ep. Parameters are n-p cbf = 0.001 and 
'^i/'^'^BB = O.OOOlSSr?^/^, where the density ratio r/ is the running independent variable. 

In Fig. 1^ we compare the mean field contributions (a) and second-order correction (b) to 
the energy functional for a ^Li, ^Li mixture, for a range of values of rj. The plots correspond 
to a situation where the scattering lengths obb = 0.2nm, aep = 2.7nm and the fermion 
density rip = 5.1 x 10^°cm~'^ are fixed, and compatible with the experiments described in 
Ref. while the boson density is varied. Note that for any reasonable boson density, the 
boson gas parameter j3 is indeed very small compared to a. In Fig. |^ we do the same for a 
^He, ^He mixture. In this case the inter-species scattering length is unknown, however we 
conjecture it to be of the same order of magnitude as the boson-boson scattering length. 
The plots correspond to a situation where a-QB = ^bf = 16nm, and rip = 3.1 x lO^'^cm"^. 
These values are compatible with current experiments on metastable triplet ^He condensates 
||T8| , p!9[l , and are particularly interesting in that the beyond mean- field corrections are quite 
large (of the order of 10%). The true significance of the boson-fermion interaction energy 
correction will of course depend on the actual value of the interspecies scattering length. We 
notice that if the latter will turn out to be of about one order of magnitude larger than in the 
pure fermionic and bosonic cases (this is for instance what happens for lithium mixtures), 
then the effect of the correction can be as large as 50% of the mean-field prediction. Then, 
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of course, also corrections proportional to the boson gas parameter have to be taken into 
account. 





FIG. 7. Boson-boson (solid line) and boson- fermion (dotted line) mean field contributions (a) 
and second-order corrections (b) to the ground-state energy-density of a '^He, ^He mixture, in 

1/3 

units of ep, to the ground-state energy-density in units of ep. Parameters are n-p cbf = 0.05 and 
n^^a-QB = O.OSr]^/^, where the density ratio r] is the running independent variable. 



C. Other physical quantities, Bose condensate depletion, and induced 

fermion-fermion interaction 



From Eq. ([70|) we can readily determine the chemical potential for the fermions /xp, 
defined as 

(dE/V\ 



/ip 



to be 



h ki 27ih aBF 

A^F = \ 

imp m 



1 H- '-^m 



The pressure reads 



5 2mp 



np + 



+ 



m 



1 + 

3n 



(75) 



(76) 



(77) 
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We then obtain for the bosonic and fermionic compressibihties, respectively, for the bosons: 



1 



dP 



47rn|aBB^^ 



and for the fermions: 



ml 



1 



27r/i^aBF 
H rip 



1 + ^-^ns) 

in 



(78) 



dP 



2 h^'ki 2-Kh^aBY 

H ub 



3 2m| 



m-pm 



(79) 



We notice that the possible instabilities induced by the mean-field boson-fermion interaction 
term in the case of a negative value of obf are contrasted by the beyond mean-field correction, 
since the latter is always positive. 

Concerning the structure of the Bose condensate fraction, besides the known depletion 
due to the boson-boson interaction, we expect in principle a further contribution to depletion 
due to the interaction of the bosons with the fermions. The depletion is computed in a 
standard way by integrating the boson propagator for the non condensed particles G'-q{p^) 
over four- momentum. To obtain the boson propagator we have to solve the Dyson equation 
(0) for G"b(p^). This yields: 



X 



P 



2mB 

Sb(p^ 



Sb(-P^. 

Sb(-p^ 



Sb(0) + Si2(0) 

1 2 



-Sb(0) + Si2(0) + 



Sb(p'^) + Sb(-p^ 



1 2 



(80) 



2mB . . 2 

where we have made use of the Hugenholtz-Pines relation. The total diagonal bosonic self- 
energy Eb(p'^) picks up a boson-boson and a boson-fermion contribution (see Eq. (|49|)). It 
can be easily checked that to first order in obb and in obf the total diagonal bosonic self- 
energy does not depend on the four-momentum. Therefore the diagonal self-energy terms 



24 



in the boson propagator Eq. ( ^0]) cancel, G'^{p^) is independent of obf, and there is no 
depletion of the Bose condensate due to the fermions to this order, since the contribution of 
the fermions to the off-diagonal self-energies vanishes anyway in ladder approximation. The 
situation will be different to next order in obf, as in this case the total diagonal bosonic 
self-energy will depend on the four-momentum. However, the calculation of Sb(p^) to 
second order in the boson-fermion scattering length at non zero four-momentum involves 
the evaluation of integrals that cannot be carried out analytically in a straightforward way. 
In conclusion, in the present situation, we will consider the depletion due to the bosons only, 



which is well known ||25| , |2^ : 



8 n-Baj , . 

riB - no = -\ ub ■ (81) 

O V TT 

We now turn to the discussion of the fermion-fermion interaction induced by the presence 
of the bosons. Subtracting the bosonic contribution from the energy density, we get 



E 



377(1+5) 37r(l + 5) TT 



(82) 



This describes the first three terms of a power expansion in a, of exactly the same form as 
that of an imperfect (unpolarized) Fermi gas, though clearly with different coefficients. There 
is thus, as expected, an induced fermion-fermion interaction, which can now be computed 
by exploiting the expressions that we have derived for both the bosonic and the fermionic 
chemical potentials. This will yield a modification of the known induced fermion-fermion 
interaction previously discussed in mean field approximation |T^. The expression for the 
induced interaction at zero energy- momentum transfer U\nd{(l^ = 0) is : 



uuq' = 0) 



dn 

which in our case reads: 



/ \ 2 



B 



(83) 



The extension to finite momentum transfer is achieved by introducing the boson density- 
density response function x{l^)i and is represented diagrammatically in Fig. where two 
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fermions interact by exchanging a boson density fluctuation wave. This is the simplest 
diagram by which two fermions can interact via the exchange of a bosonic excitation. 



FIG. 8. Effective fermion-fermion interaction due to exchange of boson density fluctuations. 

The density- density response function xil^) is independent of obf to first order. This 
can be easily verified following the same line of reasoning described previously in the analysis 
of the bosonic propagator. Thus its expression is the same as in the pure bosonic case (to 
this order): 



The vertex gpp can be determined by considering the limiting expression for xiQ^) 



xiqn 



nnoq'^/mB 



(85) 




Then the induced interaction potential in the static case (g° = 0) reads: 



f/ind(0,g) 



Aiih^al^ (1 + 4/((5)aBFA;F/3vr)^ 



(87) 



In real space this is: 



f/ind(0,g) 



(88) 




where 



1 



(89) 
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We observe that, compared to the mean-field result [0, while there are quantitative mod- 
ifications in the pre-factors, there is no qualitative change in the form of the induced in- 
teraction, i.e. we still have an attractive Yukawa potential. Modifications in the analytic 
form of the induced fermion-fermion interaction potential will appear only once second-order 
effects in the boson-fermion scattering length obf to the depletion of the Bose condensate 
are included. 



VI. DISCUSSION AND OUTLOOK 

In summary, we have determined ground-state properties of a homogeneous system of 
bosons mixed with spin-polarized fermions at zero temperature. We have calculated the 
boson-fermion T-matrix and the corresponding self-energies. Then we have shown how to 
incorporate the effects of the boson-boson interaction and derived some relevant physical 
quantities of the system, in particular the ground-state energy. The importance of the be- 
yond mean-field corrections has been discussed in several different instances of experimental 
interest. For mixtures of bosonic and fermionic Helium we have shown that the beyond 
mean-field terms may yield significant corrections (up to 50% of the mean- field result). We 
have provided partial results also on two very significant physical quantities, namely the 
Bose condensate fraction and the induced fermion-fermion interaction. To provide more 
quantitative predictions for these quantities, as well as for the BCS transition temperature, 
we will need to compute in detail the corrections to second order in the boson-fermion scat- 
tering length. Results of this analysis, which goes beyond the scope of the present paper, 
will appear in a forthcoming work, together with a detailed numerical analysis of the con- 
ditions for stability and for phase separation. Collective modes, effective fermion mass, and 
excitation spectra fully evaluated to second order in the boson-fermion scattering length will 
be discussed as well. 

Extensions of the formalism developed in this paper can be made in different directions. 
First, one can consider unpolarized spin- 1/2 fermions. Calculations are very similar to 
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the present situation with the main difference that one has to include the effects of the 
direct interactions of fermions with different spins. This would correspond to having a third 
scattering length off- As in the previous case, we expect that in the two-particle scattering 
approximation the energy contribution for a pure fermion system of spin-1/2 fermions is 
simply added to Eq. ([7^) in this case (and of course, the Fermi momentum has to be 
modified appropriately). 

The formalism can be also extended to consider finite temperature. As in the case of 
pure bosons we expect considerable difficulties near the critical temperature. Well below 
that temperature, however, we expect no major complications and the calculations will be 
similar to the present case, except that boson loops will have to be taken into account and 
frequency integrals will have to be replaced by Matsubara frequency sums. 

A very important and natural possible extension is the investigation of inhomogeneous, 
e.g. harmonically trapped, systems. This can be done by augmenting the existing mean 
field calculations via the correlations terms in local density approximation. To this end, 
the results obtained in the present work are needed. The method and the full numerical 
procedure will be described in a forthcoming paper. 

Finally, higher-order corrections may be in principle computed. These higher-order terms 
will involve also the bosonic gas parameter as well as three-particle correlations, and thus 
expansions like Eq. ([70|) will not reduce to sums of terms, where only one scattering length 
appears at a time, but will include also terms which contain products of powers of both 
scattering lengths. To even higher orders non-universal properties like the parameters de- 
scribing the shape of the interaction potentials will become important and will have to be 



taken properly into account, as it has been recently done in the pure bosonic case |31,45 
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APPENDIX A: NORMAL ORDERED PRODUCTS AND THE VACUUM STATE 

If we expand the field operators in terms of momentum eigenstates we get: 

$(x) = ^ E "ke^*^- = -^do + 0(x), (Al) 
^(x) = ^^fe^e^k-x ^ ^/(x) +4(x), (A2) 

vv k 

with 

^i^) = ^7 E «ke^^-, (A3) 

V V/ |k|>0 

V |k|<fcp 

^2(x) = ^7 E ^ke^*^- . (A5) 

In terms of the bosonic and fermionic occupation number operators /VB(k) = a^Qk and 
7VF(k) = bl.bk the ground state of the non-interacting system |^o) can be characterized by: 

Nb{0)\Co)^Nj,\Co), (A6) 

iVB(k)ieo) = Ofor |k| >0, (A7) 

^F(k)|^o) = |eo)for |k| <A;f, (A8) 

iVF(k)lCo) = Ofor|k| > fcp, (A9) 

where A^b is the total number of bosons, which in this case coincides with the number of zero- 
momentum bosons A^o (Bose- Einstein condensate). In occupation number representation we 
thus have 
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I^o) = lA^B, 0, 0, . . . >_B 1, . . . , 1, 1, 0, 0, . . . >F , (AlO) 

where the subscript B refers to the boson Hilbert space and the F to the fermion Hilbert 
space. The change from 1 to in the fermion state happens at kp. Additionally, 

0(x)|eo) = ^i(x)leo) = ^2(x)|^o) = 0. (All) 

In this sense the ground state can be regarded as the vacuum state with respect to fermions 
excited above the Fermi sea, the fermion holes below the Fermi sea, and the non-condensate 
bosons. 

The normal product is defined on pairs of creation and destruction operators: 

: V^,(x,t)ViI(x',t') : = -V^I(x',t')^.(x,t), 
:0(x,t)^j(x',t'): =^J(x',O0(x,t), 

: ^,(x,t)0t(x',t') : = 0t(x',t')^^.(x,t), (A12) 

for j,k G {1,2}. For all other pairs of creation and destruction operators the normal product 
is the same as the ordinary operator product. It can also be readily determined that 

[0(x,t),0t(-/^t')] = (^o|0(x,t)0t(x',t')|^o), 

{^i(x',t'),v^i(x,t)} = (eoi^Hx,t')^(x,t)ieo), 

{^Pix,t),ij^ix',t')} = (eo|^(x,t)^t(x',t')|^o), (A13) 



and all other (anti-) commutators are zero. With Eqs. [A12| and |A13| , the contractions of Eq. 
(P7|) can be readily evaluated. 

APPENDIX B: EVALUATION OF THE T-MATRIX AND COUPLING 
CONSTANT RENORMALIZATION 

1. The first integral I 

We define 
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^ ^ / g(|P/2-k|-fcp) 

J nP^-n\P/2 + k)y2mB-h\P/2-k)y2mF + H + iiy' ^ ' 

Transforming the integration variables to P/2 — k gives: 

J h\^/2m-h^P -k/mB-hPO + h^P^/2mB- fi-ii^ 
Setting a = h'^/2m, b = fi^P/m-Q and E = —hP'^ + h'^P'^ /2'm-Q — fx and transforming to 
spherical coordinates we get: 

X = 2tx i dkk"^ / dd) sin 



F 







ak"^ — hk cos 6 + E — 



= — / dfcfcln— , B3 

b Jkp ak^ + bk + E-w ^ ^ 

where we will ultimately consider the limit k^ oo. Using D = (fe/2a)^ — E/a we can 

approximate for small u {if D ^ 0; the case D = can be treated similarly and gives the 

same answer as taking the limit D ^ at the very end): 

h^p 



X In 



F 

k + mP/niB + VD + iu/2aVD 
k — mP/m-Q — \fD — iv /2a\/15 



, k + mP/rriB- \fD -ip/2a\fD\ 

+ m 1= !=]■ (-d4) 

k - mP/mB + VD + iu/2aVD J 



The integral can be solved ^J] to give 

STimkc AirrnkF 
lim X = 2 1 2 — 



TT (rriBkl m^P r— mBD\ 

+ 72 2mVD — 

\ P niB P J 

kp + mP/niB + VD + iu/2a\^ 



X In 



kp — mP/rriB — — ii'/2a^/D 

71 fniBkl rn^P r— mBD\ 
-0 — 7^ 1- 2my D -— 



, k^ - mP/niB + \fD + iv l2a\fD 

X In ■= -=, (B5) 

kp + niP/niB — y D — iv /2ay D 

where outside the logarithms we have taken the limit z/ ^ (simply setting i/ = 0), and we 

have made use of the identity 
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2i 1 + 
lim X 111 



2ax, 



(B6) 



1 — a/x 

for the limit kc ^ oo. There remains an ultraviolet divergent term; the boson-fermion 
T-matrix [Eq. (0)] is however ultimately renormalized by the second integral. 

The real part of X is readily evaluated in the limit z/ ^ by setting i/ = and using the 
absolute values inside the logarithms: 



lim ReX 



+ 



71 f m^kp rn^P m^D\ 



X In 



(A;f + mP/rriBf - D 



27rm 



Din 



{kp + vHjf - {mP/m-Bf 



(B7) 



Using the identity (easily evaluated by polar decomposition) 

a + iu : sgn(a) = sgn(6) 
lim Im In = < , 

TT : sgn(a) ^ sgn(6) 
the imaginary part of X in the limit z/ — > can be evaluated to be: 

lim ImX = ^ f - !^ - 2mv^ - 

i^-O \ P TUB P I 

if > and fcp < \mP/mB - v^|; 



lim ImX 



2 ' 



if D > and \mP /tub — \fP>\ < kp < mP/niB + ^/D', and 



(B8) 



(B9) 



(BIO) 



lim ImX = 0, 



(Bll) 



if D < or fcp > mP/niB + VD- 
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2. The second integral J 



We define 



a Jo \ K'^ — kf — If /a ^ 
where as before a = /2m, and we liave transformed to polar coordinates and integrated 



-I dk l-—^-^^), (B12) 



over tlie angle variables. The integral can be evaluated MM to give 



Ank^ 2n r- -— kc + ^kf-iu/a 

J = \ kf — IP I am V = 

« « K-^kl-iv/a 



We then use 



to get 



27r I x/ki—iu/a 

+ fl B-tu/a\n ^ / =. (B13) 



lim In ^ = (B14) 

kc — Jkf — iu/a 



limj=— B15 



If we now take the sum of Eqs. (p5|) and ( pi5|) , the ultraviolet divergent terms cancel exactly. 



The resulting expression for X + ^7 can then be substituted into Eq. (^) to get Eq. (|59|) for 
the renormalized boson-fermion T-matrix. 
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